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main rotor and control surface lift-curve slope, 1/rad 

main rotor disk area, ft 2 

main rotor and control surface chord, ft 

section drag coefficient 

lift coefficient 

coefficient of thrust in ground effect 

coefficient of thrust out of ground effect 

c.g. location from base of pilot compartment, ft 

vertical distance from base of helicopter to hub, ft 

vertical distance from base of helicopter to pilot c.g., ft 

vertical distance from hub to rotor tip, ft 

diameter of propellers, ft 

total drag, lb 

profile drag of the main rotors, lb 
profile drag of the control surfaces, lb 
differential perturbation gain, 1/sec 2 
differential perturbation gain, 1/sec 2 
hinge offset (% main rotor radius) 

control surface lift force of the rotor corresponding to \|/|=o = 0°, lb 

control surface lift force of the rotor corresponding to \j/ t =0 = 1 80°, lb 

gravitational acceleration, ft/sec 2 

height from ground to average vj location, ft 

rotor tip height of the rotor corresponding to \j/ t= o = 0°, lb 
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hi rotor tip height of the rotor corresponding to \|/ t= o = 180°, lb 

h 0 initial height of the base of the helicopter, ft 

hR distance between c.g. and hub, ft 

lb mass moment-of-inertia of one main rotor, slugs-ft 2 

lx roll axis mass moment-of-inertia, slugs-ft 2 

I y pitch axis mass moment-of-inertia, slugs-ft 2 

K a differential height-to-angle gain of the control surface corresponding to h r rad/ft 

Kb differential height-to-angle gain of the control surface corresponding to hi, rad/ft 

K Cj i roll feedback gain to roll acceleration, 1/sec 2 

K c ,2 pitch feedback gain to roll acceleration, 1/sec 2 

Kd,i roll feedback gain to pitch acceleration, 1/sec 2 

Kd,2 pitch feedback gain to pitch acceleration, 1/sec 2 

l C s length of control surface, ft 

L v lateral velocity stability, rad/sec/ft 

Lp roll damping, 1/rad/sec 

m total mass of the helicopter, slugs 

mb mass of the pilot and pilot compartment, slugs 

m r mass of one main rotor, slugs 

m t mass at one rotor tip, slugs 

Mq pitch damping, 1/rad/sec 

M u longitudinal velocity stability, rad/sec/ft 

N number of blades 

N r yaw damping, 1/rad/sec 

N v directional velocity stability, rad/sec/ft 

p body axis roll rate, rad/sec 
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roll acceleration feedback, rad/sec 2 

body axis pitch rate, rad/sec 

pitch acceleration feedback, rad/sec 2 

body axis yaw rate, rad/sec 

radius of main rotors, ft 

effective radius of main rotors, ft 

contol surface area of one control surface, ft 2 

thickness of main rotors at 30% chord, ft 

thrust required in ground effect, lb 

thrust required out of ground effect, lb 

time, sec 

body axis longitudinal velocity, ft/sec 
body axis longitudinal wind gust, ft/sec 
body axis lateral velocity, ft/sec 
induced velocity in ground effect, ft/sec 
induced velocity out of ground effect, ft/sec 
body axis vertical velocity, ft/sec 
body axis vertical wind gust, ft/sec 
weight, lb 

change in X-force with respect to q, ft/sec/rad 

longitudinal damping, 1/sec 

change in Y-force with respect to p, ft/sec/rad 

lateral damping, 1/sec 

vertical damping, 1/sec 

beta, coning angle, rad 


Sr.cs rate and position limited control surface deflection of the rotor corresponding to \j/t=0 = 0°, rad 

Si,cs rate and position limited control surface deflection of the rotor corresponding to \i/t=0 = 1 80° 

rad 

8 r.cs control surface deflection signal to actuators of the rotor corresponding to y t= Q = 0°, rad 

8 i >cs control surface deflection signal to actuators of the rotor corresponding to y t =0 = 180°, rad 

8cs,d differential control surface deflection, rad 
y gamma, Lock Number 

£2 omega, rotor speed, rad/sec 

<t> phi, body axis roll angle, rad 

¥ psi, body axis yaw angle, rad 

p rho, air density at sea level, slugs/ft 3 

o sigma, solidity 

0 theta, body axis pitch angle, rad 

00 theta zero, main rotor initial pitch angle, rad 

X inflow 
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SUMMARY 


This report documents the study of a control system for the Da Vinci II human-powered helicopter 
in hovering flight. This helicopter has two very large, slowly rotating rotor blades and is considered to be 
unstable in hover. The control system is designed to introduce stability in hover by maintaining level 
rotors through the use of rotor tip mounted control surfaces. A five degree of freedom kinematic model 
was developed to study this control system and is documented in this report. Results of this study show 
the unaugmented configuration to be unstable due to the large Lock Number, and the augmented configu- 
ration to be stable. 

The reason for NASA’s involvement in this study (and the publication of this document) was so 
that instructors and students at the university level would have an educational aid for modeling and coding 
dynamic systems. The role of NASA in this study included the development and analysis of the kinematic 
model and control laws. Both analytical and numerical techniques were used. 


INTRODUCTION 


Since 1981 the California Polytechnic State University (Cal Poly, San Luis Obispo) student chap- 
ter of the American Helicopter Society has been involved in an effort to win the Igor Sikorsky Human- 
Powered Helicopter Design Competition prize. The requirements are to achieve human-powered hovering 
flight for 1 min, to reach an altitude of 3 m, and to stay within an area of 10 by 10 m. The first prototype, 
the Da Vinci I, was built of advanced composite materials and had two 50-ft-radius rotor blades which 
tapered from an 8-ft chord at the root to a 6-ft chord at the tip. The main rotors were driven by tip- 
mounted propellers that were 6-ft in diameter and turned at 350 rpm. The pilot supplied power to the pro- 
pellers by winding-up string that was threaded through the main spars and wrapped around the shaft of the 
propellers. Rotor speeds up to 6 rpm could be obtained by this prototype. 

The Da Vinci II differs from the Da Vinci I in that it features two 67-ft rotor blades having constant 
chords of 3 ft, refined advanced composite technology, tension cable reinforcement to reduce bending, and 
a unique control system concept. 

Although the Da Vinci II was designed to sustain hover for 1 min, initial flight tests of the unaug- 
mented configuration showed unstable dynamic behavior. One rotor tended to generate more lift than the 
other. The rotor generating less lift would eventually impact the ground in roughly 30 sec. The aug- 
mented configuration described in this report has not yet been flight-tested. 


AIRCRAFT DESCRIPTION 


The Da Vinci II is depicted in figure 1 . The main spars are made from carbon-graphite, filament- 
wound composites. The rotor ribs are a sandwich type construction consisting of a styrofoam core 
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covered with S-glass or graphite. The rotors are covered with Tedlar and the propellers mounted at the 
rotor ups were made by covering expandable foam with Kevlar cloth. The airfoil design of the main rotors 
is a Lissaman 7769 and the rotors are at a fixed incidence of 10°. 

The control system consists of control surfaces mounted outboard of the tip mounted propellers. 
Optical sensors are mounted near each control surface in order to measure height from the ground. These 
control surfaces are differentially driven in proportion to the difference in height measured by the optical 
sensors. The control surfaces have the same airfoil shape as the main rotors and are actuated by servos 
mounted inside the spars. 

The Da Vinci II has no tail rotor or any other conventional control mechanisms. The pilot com- 
partment is rigidly attached to the shaft and hub, as are the main rotors. 


MATH MODELS 


A block diagram representation of the kinematic model and control system of the Da Vinci II is 
shown in figure 2. A description of the axis systems is given in Appendix A. The equations and princi- 
ple assumptions used to describe both the kinematic model and the control system are presented in the next 
two sections. Special considerations and developments pertaining to the kinematic model are given in 
Appendix B. ° 


Kinematic Model 


This section contains modified, linearized perturbation equations of motion of a helicopter used as 
the kinematic model for the Da Vinci II. These equations were derived from the general equations of 
motion based on the following assumptions (ref. 1 ): 

1. The flight condition is hover. 

2. The rotors have a rectangular planform with no twist. 

3. There are no stall or compressibility effects. 

4. There is no higher harmonic rotor blade flapping. 

5. There is no pitch-flap coupling. 

6. The quasi-steady assumption is employed. 

7. The vertical, longitudinal, and lateral axes are decoupled. 

8. Small angle approximations are used. 

A detailed discussion concerning the effects of aeroelasticity, rotor tip losses, and hovering in 
ground effect on the equations of motion for the Da Vinci II is given in Appendix B. The modified, 
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linearized perturbation equations of motion based on the assumptions and special considerations are as 
follows: 
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The stability derivatives in these equations are computed from approximate relationships using air- 
craft physical parameters listed in Appendix C, table 3. Many of these stability derivatives are a function 
of I x or Iv which are a function of rotor position in the tip-path plane, \\f. This is because the large con- 
tribution of the main rotors to these inertia terms is not a constant value. The development of I x and I y is 
given in Appendix B. The stability derivative values are listed in Appendix C, table 4, and the approximate 
relationships are as follows (ref. 1): 


X u = -[T(daiNF/du) + (dHxpp/du)]/m 


where 


daiNF/du = aiNF/u = (80o/3) + 2X 
dHjpp/du = poAQRcdo/4 
X = -vi/(QR) 
vi = (see Appendix B) 

Ct = (see Appendix B) 

A = JtR 2 
c = NC/(7tR) 

T = (see Appendix B) 


X q = 
where 


-[T(daiNp/dq) + dHTPp/dq]/m 

daiNF/dq = - 16 /(yQ) 

y = paCR 4 /I b 

I b = m r R 2 /3 (inertia of a thin rod) 
dHjpp/dq = -paAo(QR) 2 ^/(2YQ) 


Y v = X u 
Y p = -X q 

Z w = -pAQR(dCx/dw)/m 
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where 


dCx/dw = aa/[8 + aaV (2 /Ct)/2] 

Ly = *M u Iy/I x 

W ^ ere lx = (see Appendix B) 

Iy = (see Appendix B) 

Lp = Mqly/I x 

= [hR(dHxpp/du + T(daiNp/du)) + dMs/du]/I y 

where dMg/du = NeR(C.F.)daiNp/du/2 

centrifugal force, lb C.F. = m r RQ2/ 2 + mi RQ2 

M q = [hR(dHxpp/dq + T(daiNp/dq)) + dMs/dq]/Iy 

where dMs/dq = -8NeR(C.F.)/(yQ) 

N v = 0 (due to lack of directional control mechanism) 

N r = 0 (due to lack of directional control mechanism) 


Control System Model 


The control system of the Da Vinci II is mathematically described in this section and all assump- 
tions and restrictions are discussed. ump 


inscription 

senceH H^r;r t0r CO T? 1 suifaces of the Da Vinci n move differentially proportional to the optically 
h f ht the mam rotor tips ' 7116 con *ol surface actuators are driven by height differ^ 

ence signals, such that these actuators increase the angle of attack of the control surface corresponding to 

senso^Thk 102 J ensor 311(1 decrease the angle of attack of the control surface corresponding to the higher 
sensor. This creates a moment proportional to the measured height difference. The optically sensed § 

height and height-difference signals transmitted to the control surface actuators are considered to be accu- 

each^t 111 ' TSmrJ thC PUrPOSCS ° f tHiS StUdy ' The C ° ntr01 SUlfaCeS move linearl y one d ^e for 

each foot of height difference measured (this is discussed further in the Numerical Method section) 


Mathematical Development 


The mathematical development of this control system utilizes small-angle 
consists of three events described by equations (1) through (10): 


approximations. It 
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1 . Determination of optically-sensed height, height difference, and actuator signals. 
Optically Sensed Height: 

h r = h Q + d a + d c - wt - R[sin $ cos Y * sin 9 sin yl 

hi = h Q + d a + d c - wt - R[sin $ cos(\|/ + it) - sin 0 sin (y+ it)] 

Height difference: 

h r - hi = -2R[<J) cos Y - 0 sin \j/] 

Actuator signals: 

8 r iC s = (hr " hi)K a 

8'l,cs = (hr * hi)Kb 


( 1 ) 

( 2 ) 

(3) 

(4) 

(5) 


. v = K 1 (deg/ft) 

where Kb K a 57.3(deg/rad) 

2 . Calculation of the rotational accelerations about the c.g. generated as a function of the resultant 
control surface lift forces. 

Resultant control surface lift forces: 


F r = -Fi = -0.5p[O(R + lcs/2)]^S cs a 8 r ,cs 

where 8 r , C s = rate and position limited value of 5' r>C s 

The analytical method presented in this report evaluates the control system using 5' r ,cs, 
whereas the numerical method uses 5 r)CS . The implications of this will be discussed in the 

next section. 

Generated rotational accelerations: 


pfb = 2F r (R + W2)cos y/Ix 
= K c ,i<|) + K c ,20 

where K c i = 2RaQ 2 pK a (cos t}0^Scs(R + lcs/2) 3 /Ix 

K Ci 2 = -2RaQ 2 pK a cos y sin \| / S cs (R + W 2 ) 3 /Ix 

qfb = -2Fr(R + lcs/2)sin 
= Kd,i<() + Kd,20 
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where 


Kd,l = -2Raft 2 pK a cos \|/ sin \| / S cs (R + l cs /2) 3 /I y 
Kd,2 = 2Raf2 2 pK a (sin y) 2 S cs (R + l cs /2) 3 /I y 
3 . Feedback of the generated accelerations to the helicopter body-axis accelerations. 

p = L v v + Lpp + KcjtJ) + K C 20 ^ 

q = M u u -I- Mqq + K di i<f> + K d 2 0 (1 

contr °l system is designed to drive the measured height difference to zero by generating a 
restoring moment and associated acceleration. The restoring moment is only present when height differ- 
ences are present. Contributions of the control surfaces to thrust and induced velocity are neglected and 
are discussed in Appendix B. 6 


ANALYSES 


. testability characteristics of the Da Vinci II have been studied using analytical and numerical 
methods. The analytical method entailed development of root locus plots in order to define sources of 
instability of the unaugmented configuration as well as to determine the effect of the control system on the 
stability of the Da Vinci II in hover. The numerical method entailed development of a discrete simulation 
tor use as a design tool suitable for determining appropriate control system design specifications (e.g. 
actuator rate limit, actuator position limit, and control surface area). 

The kinematic model previously described is decoupled in the longitudinal, lateral, and vertical 
axes. However, the introduction of the control system couples the longitudinal and lateral axes. This is 

rntnr?^ 6 3Ct “ at10 " of t ? e c° ntr °l surfaces can induce accelerations in both roll and pitch when the main 
rotors are not aligned with the x or y axes. 

Wind gust perturbations were used in order to study the response of the unaugmented and aug- 

^l n Hl h T flgUra r, S Under similar conditions. The block diagram given in figure 2 depicts the location 
at which these perturbations are introduced to the kinematic and control system models. A differential 
control-surface deflection input has also been depicted in figure 2 as an alternative perturbation to the 
kinematic and control system models, but this perturbation was not used in this study. It has been 
depicted for illustrative purposes only. The perturbation gains associated with the differential control 
surface deflection input are defined by equations (11) and (12). 


Dq - p£2 2 (R + l C s/2) 3 S cs a sin \|//I y 


( 11 ) 


D<j) - -pf2 2 (R + l cs /2) 3 S cs a cos \j//I x 
The resultant mathematical equations and perturbations are as follows: 


( 12 ) 
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where u = X u (u + u') - g9 + Xq0 


0 = M u (u + u') + + Mq0 + Kd,i<t> + De5 cs ,d 


Analytical Method 

Root locus plots were developed based on the kinematics and control system equations describing 
the Da Vinci II. The goals were to define any instabilities of the unaugmented configuration and study the 
effects of the control system on the stability of the Da Vinci II in hover. 

Characteristic polynomials for the unaugmented and augmented configurations were developed and 
are based on the use of perturbations as inputs and aircraft states as outputs. The development of the 
characteristic polynomials is given in Appendix D and associated root locus plots are depicted in figures 3 
and 4 (ref. 2) for values of control surface area obtained from the results of the discrete simulation, dis- 
cussed in the next section. 

The plot depicting the augmented configuration was developed using 5' r , C s ( not 8r iCS ) in the con- 
trol system feedback equations previously described. This is because rate and position limits represent 
nonlinearities in the modeled system which cannot be meaningfully represented by root locus analysis. 
Therefore, the plot depicted in figure 4 is not truly representative of the actual control system model but 
does illustrate the stability characteristics of the control system with no rate or position limitations. 

The control system has no effect on the stability of the vertical axis and it remains decoupled from 
the longitudinal and lateral axes. The vertical axis is described by a first-order-lag, and is stable because 
Z w is negative, as given by equation (14): 

w(t) = Z w w(t) + Z w w'(t) 

w(s) z w (i4) 

w'(s) s - Z w 

The longitudinal and lateral axes are decoupled for the unaugmented configuration, and are 
described by the characteristic polynomials given by equations (15) and (16): 
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(15) 


1 - M,, — ^ qS ' 

u s(s - X u )(s - M q ) 


= 0 


(16) 


IT ( Y P S + fi) n 
Lv s(s- Y v )(s-L p )-° 

TTie longitudinal and lateral axes are coupled for the augmented configuration, and are described by the 
characteristic polynomial given by equation (17): * 

1 K c . 2 K d ,i(s - X u )(s - Y v ) 

{ (s " x u)[s(s - Mq) - Kd,2l - M u (X q s - g)J {(s - Y v )ls(s - L p ) - K Cjl ] - L v (Y p s + g) ) “ 0 

(17) 

The stability derivatives and control system gains M q , M u , Lp, L v , K^, K^, K dt i, and K d 2 are func- 
tions of V|/. Thus, the stability characteristics of the longitudinal and lateral axes vary as a function of y, 

3.S WCIl, 


Root locus plots of the unaugmented configuration of the Da Vinci II are depicted in figure 3 for 
the longitudinal and lateral axes for values of y = 0°, 45°, and 90°. Roots corresponding to all other val- 
ues of v vary between the roots at y = 0° and V = 90°. Regardless, for any given value of w two 
complex poles and one real zero are in the unstable region. The unstable poles correspond to s and 
(s - Mq) for the longitudinal axis, and s and (s - Lp) for the lateral axis. The exact location of the poles 
and zeros depicted by the root locus plots in figure 3 are given in table 1. 


TABLE 1.- POLES AND ZEROS OF UNAUGMENTED CONFIGURATION 


Longitudinal Axis Lateral Axis 

\|/ = 0° 

zero 

Poles 

18.94 

0.078 ±1.4i, -12.08 
0.73 

Zero 

Poles 

Ly 

18.94 

0.014 + 0.59i, -11.86 
-0.13 

1 V = 45° 1 

Zero 

Poles 

M u 

18.94 

0.027 ±0.79i, -11.97 
0.23 

Zero 

Poles 

L v 

18.94 

0.027 ±0.79i, -11.97 
-0.23 

L \j/ = 90° 1 

Zero 

Poles 

Mu 

: 18.94 

0.014 ±0.59i, -11.86 
0.13 

Zero 

Poles 

Ly 

18.94 

0.078 ±1.4i, -12.08 
-0.73 


The reason two complex poles are in the unstable region for the longitudinal and lateral axes is 
because the values of M q and Lp are very small. This places the poles (s - M q ) and (s - Lp) very close to 
he origin i of each root locus plot, next to a pole at the origin, for values of M u and Lv equal to zero. The 
locus of these poles quickly diverge toward the zero in the unstable region as M u and L v vary to their 
actual values. The common denominator in the M q and Lp terms which make them so small is the Lock 
Number, y (ratio of aerodynamic forces to inertial forces), which is the ultimate cause of the instability of 
the unaugmented configuration for the given values of M u and L v . 
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A root locus plot of the augmented configuration of the Da Vinci II is depicted in figure 4 for va - 
ues of y = 0°, 45°, and 90°. Roots corresponding to all values of y vary between the roots at y - 0 
and \i f = 45°. Regardless, for any given value of V at least two real poles are in the stable region. 
Additionally, one complex pole pair varies between the unstable and stable regions, becoming stable as V| / 
approaches 45°, again at 135°, 225°, and so on. Another complex pair remains in the unstable region and 
moves away from the origin as V approaches 45°, again at 135°, 225°, and so on. The exact location o 
the poles and zeros depicted by the root locus plot in figure 4 are given in table 2. 


TABLE 2.- POLES AND ZEROS OF 
AUGMENTED CONFIGURATION 


< 

II 

o 

o 

Zeros 
Poles 
Poles 
Poles 
K c .2 Kd.l 

-11.81, -11.81 
-11.81, -12.12 
0.014 ±2.32i 
0.08 ±1.40i 

0 

II 

o 

Zeros 

Poles 

Poles 

Poles 

K c .2 Kd.l 

-11.81, -11.81 
-11.73, -12.05 
1.105 ±2.79i 
-1.059 ±2.79i 
36.13 

< 

ii 

s 

o 

Zeros 
Poles 
Poles 
Poles 
K c .2 Kd.l 

-11.81, -11.81 
-11.81, -12.12 
0.014 ±2.32i 
0.08 ±1.40i 

0 


A close inspection of the root locus plot depicted in figure 4 reveals that at least one of the unstable 
complex pole pairs remains close (less than or equal to 0.014) to the imaginary axis in the unstable region 
for \i f = 0°, 90°, 180°, and so on. The complex poles closest to the imaginary axis are considered to be 
dominant if the ratio of the real parts of these poles to the real parts of the next closest poles are greater 
than five (ref. 5). This is the case for the augmented configuration at y = 0°, 90°, 180°, and so on 
because the real parts of the unstable complex poles have a ratio of approximately six to the next closest 
roots. However, as y approaches 45°, 135°, 225°, and so on the dominant unstable complex poles move 
away from the origin and lose their dominance due to the presence of complex poles in virtually the same 
location in the stable region. The augmented configuration essentially becomes stable at these values of y 
because of the cancellation effect of the unstable complex poles and the presence of the stable, real poles. 
The characteristics of the augmented configuration vary from slightly unstable to stable for different values 
of y. This analysis represents the characteristics of the control system with no rate or position limiting. 


Numerical Method 

A discrete simulation was developed to model the kinematics and control system of the Da Vinci II. 
A description of this simulation is given in Appendix E. Results of the simulation are given in the form of 
time histories of the state variables in the presence of the wind gust perturbations described previously. 
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control Ilrl 7 Stei J senses P r Siti ° n and produces accelerations (pn> and q ft ) through the use of the 
control surfaces. Rate and position limits are physical constraints of the actuator mechanisms used to 

c ange the angle of attack of the control surfaces. The rate and position limits were studied in the discrete 
simulation because they produced lags and decreased the authority, respectively, of the control surfaces. 

e goal was to reduce the effectiveness of the control surfaces in responding to differential rotor height so 
t at the accelerations produced would not be extreme thereby creating an unstable system. Corre- 
spondmgly, caution was used in sizing the limits so that the control surfaces would not produce 
insufficient accelerations. Control surface area was also varied in the discrete simulation in order to 
configur^o * ^ anc * thereby changing the stability characteristics of the augmented 

Various combinations of rate limit, position limit, and control surface area values were studied for 
the control system with a step input of a 5-mph forward velocity wind gust. This was a reasonable distur- 
bance because the intent was simply to excite the augmented configuration with the same perturbation used 
to induce large, unstable motions from the unaugmented configuration. A lateral-velocity wind gust per- 

conditicm Ys^over ^ pr ° duCed the Same results because the configuration is symmetrical and the flight 

J° quantitatively indicate the effectiveness of the control surface parameter being evaluated the roll 

Hi’viHpHh C abs ° l fl Va 7f S T re summed ever y cycle through the duration of the perturbation and 
ded by the sum of the roll and pitch angle absolute values ten seconds later (one rotor revolution after 

aS “gl 5°Lc 5 ^ dUra “° n - ThC ra,i ° “ defmed be ' OW fOT 2 SKP inpU ‘ introduced 


I(l<t)l + 191) 

2(i<|)i + iei) 


(from t = 1.25 to 3.75 sec) 
(from t = 11.25 to 13.75 sec) 


(18) 


1 Vidues s reaterthan one indicate stability, values equal to one indicate marginal stability and values 
oTLrt^tTri indlcat ® instablll [y because, respectively, progressively smaller motion changes as the result 
of perturbations would make the denominator of the above equation smaller than the numerator no differ- 

mofion S° n Chan w S WO l uld res . ult ln the denominator equalling the numerator; and increasingly larger 
motion changes would make the denominator larger than the numerator. The process used to obtain 

S‘ a 7 rate i™ 1 ’ positlon and contro1 surface area values was to vary one parameter while the 

f?h Sta 7 ^ Va UC ° f thC p?rameter yieldin g the hi g h est ratio was chosen and held 
constant while one of the other two parameters was varied. This iterative process was repeated several 

times until a maximum ratio value was obtained. Plots for this process are given in figure 5 for parameter 
vrnanon about values obtained from the final iteration. The final value for S* of on! control surface is 
115 ft the final value for the rate limit (R.L.) of each control surface actuator is 0.24 rad/sec, and the 
final value of the position limit (P.L.) of each control surface is 0.04 rad. 

at 0 Ol^trad/fr^ ° f ** ?* back gai "’ K * WaS ch osen and held constant 

toii 745 V dg f f) L f ° r he purposes of thls study. This parameter represents the amount of con- 
trol surface angle of attack obtained per foot of rotor-tip-height difference. Varying K a is mathematically 

equivalent to .varying the control surface area, as can be seen from the controlsystem equations presented 
previously. The same responses would have been achieved for different values of S^and K a as long as 
the product (R + lcs/2) 3 Sc S K a remained equal to -6.57 x 10 4 ft^-rad (where l cs = S cs /C). The value of § this 
product was used in generating the root locus plot described in the previous section. 
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Unaugmented time histories are given in figure 6 and augmented time histories without rate and 
position limiting are given in figure 7. These numerical results agree with the analytical results presented 
previously. The root locus plots given in figures 3 and 4, and the time histories given in figures 6 and 7, 
show that the augmented configuration with no rate or position limiting improves the stability charac- 
teristics of the Da Vinci II when compared to the unaugmented configuration. Exact correlation between 
the analytical results and numerical results cannot be determined because the system is nonlinear and the 
natural frequency and damping of each configuration varies as a function of \j /. Additionally, time histo- 
ries of the augmented configuration with rate limit and position limit values of 0.24 rad/sec and 0.04 rad, 
respectively, are given in figure 8. 

Perturbations are introduced at \j/ = 45° (t = 1.25 s) and last until \y = 135° (t = 3.75 s). The 
results depicted in figure 6 show oscillatory instability of the unaugmented configuration. The positions, 
rates, and accelerations depicted in figures 7 and 8 are substantially less than those depicted in figure 6 and 
eventually subside to zero. The positions and rates depicted in figure 8 appear fairly well damped, except 
for the roll rate response, when compared to the positions and rates depicted in figures 6 or 7. However, 
the maximum roll rate attained never exceeds 10 deg/s and subsides in roughly ten seconds as well. The 
accelerations depicted in figure 8 do not appear well damped. This behavior is due to the summation of the 
rotational accelerations produced by the control surfaces (pn> and qfb) with the body axis roational 
accelerations whenever a differential rotor height is present. The control surfaces are extremely active in 
the first ten seconds and are sensitive to rotor tip height difference (h r - hi) even with rate and position 
limiting. However, the maximum pitch and roll accelerations obtained are no greater than 12 deg/sec 2 and 
0.8 ft/sec 2 , respectively. 


RESULTS 


Three conclusions can be drawn from the results of the analyses of the kinematic model and control 
system, presented in this report. The first conclusion is that the unaugmented configuration is unstable, 
and this is evidenced by the root locus plots depicted in figure 3 and the time histories depicted in figure 6. 
The second conclusion is that the control system without rate and position limiting improves the stability 
characteristics of the Da Vinci II, and is considered slighdy stable when compared to the unaugmented 
configuration, and this is evidenced by the root locus plot depicted in figure 4 and the time histories 
depicted in figure 7. The third conclusion is that the control system with the stated combination of values 
of the control system design parameters (K a , S cs , lc S , rate limit, and position limit) is a stable system. This 
is evidenced by the time histories depicted in figure 8. 

The Da Vinci II is a second generation human-powered helicopter prototype. There is no actual 
flight test data available to validate the kinematic model which was used to study the control system. 
Although the results depicted in figure 8 show the control system is able to stabilize the Da Vinci II in 
hover, these results must be considered preliminary until the kinematic model is validated by comparison 
with actual flight test data or until proven by flight tests. 
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APPENDIX A 


AXIS SYSTEMS 


The equations describing the kinematics are used to calculate the accelerations about the c.g. of the 
Da Vinci n. These are body axis accelerations. The body axis system has its origin about the c.g. and is 
depicted in figure 9. 

The equations describing the control system use ground height to calculate control surface deflec- 
tion. The ground height is referenced relative to the earth. The earth axis reference system used is also 
depicted in figure 9. 
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APPENDIX B 


SPECIAL CONSIDERATIONS 


This appendix contains the analysis and development of mathematical relationships which describe 
elements of this study which are unique. There are four elements and associated assumptions, as follows: 

L Aeroelasticity: affects the calculation of thrust and mass moment-of-inertias. 

2. Tip-losses: reduce the effective rotor radius, however the effect is negligible. 

3. Ground effect: affects the calculation of thrust and induced velocity. 

4. Control surface contributions to thrust and induced velocity: are considered negligible. 

The justification and associated mathematical development for these assumptions are herein described. 

The first consideration is aeroelasticity. The Da Vinci II was designed with tension cables to 
reduce bending of the main rotors in hover. The measured tip deflection and coning angle, P, during a 
flight test were roughly 12 ft andlO°, respectively, and no main rotor pitching or flapping motion was 
noted. Based on this information the Da Vinci II was modelled using rigid-rotor equations at a constant 
coning angle of 10°. The inertial representation is depicted in figure 10, and the calculation of the c.g. 
location, I x , and I y are given by equations (19), (20), and (21), as follows: 

1 . Center-of-gravity calculation: 

m d = 2[m r (d a + dc/2) + mt(d a + d c )] + m^d^ (19) 

d = 5.5 feet 

2. I x calculation: 

I x = 2mbdb 2 /3 + mb(d - db) 2 

(contribution of pilot and pilot compartment) 

+ 2mr(R cos P cos \)/) 2 /3 + 2m r (d a + djl - d) 2 
+ 2m r (R cos P sin y) 2 /12 
(contribution of main rotors) 

+ 2m t (R cos P cos \j/) 2 + 2m t (d a + d c - d) 2 
(contribution of mass at rotor tips) 

I x = 221.0 + 1133.3 (sin \|/) 2 + 7253.2 (cos \j/) 2 (20) 
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3 . I y calculation: 

I y = 2mbdb 2 /3 + mb(d - db) 2 

(contribution of pilot and pilot compartment) 

+ 2m r (R cos p sin \|/) 2 /3 + 2m r (d a + dc/2 - d) 2 
+ 2m r (R cos p cos y) 2 /12 
(contribution of main rotors) 

+ 2m t (R cos p sin \|/) 2 + 2mt (d a + d c - d) 2 
(contribution of mass at rotor tips) 

I y = 221.0 + 1133.3(cosi|/) 2 + 7253.2(sin\|/) 2 (21) 

Thrust is also affected by aeroelasticity because the resultant lift vector of each rotor is tilted inboard by the 
amount of the coning angle. The thrust required to hover is calculated by equation (22), as follows: 


Too = 



W 2 + Dp ' 2 
(COS P ) 2 


( 22 ) 


where D 0 = CdopCR(£2R) 2 - 8.5 lb 

W = 285 lb 
Too = 290.47 lb 

The second consideration is tip losses, which tend to reduce the effective rotor radius. Equa- 
tion (23) can be used to define the effective rotor radius, as follows (ref. 3): 


Re, 

R “ 1 " N 


(23) 


where 


p iOQ 

To ° “ pA(QR) 2 
Ct«= 0.00489 


For this study equation (23) is modified to account for coning which reduces Rg by cos P, and is given 
by equation (24), as follows: 


Re = ' 


>/ CTo. 


N 


Rrcos p 


(24) 


R e = {64.66}cos p = 63.86 ft 

The effective radius, Re, calculated above constitutes a 4.7% reduction in size of the actual radius, R. 
Ground effect is estimated to decrease the value of the thrust coefficient by 16.7% as is justified further on 
in this appendix. This decrease in the thrust coefficient increases Rg to 64.06 ft. Furthermore, the 6-ft 
diameter, tip-mounted propellers generate an induced velocity component perpendicular to the rotor span, 
which would tend to increase Re to an even greater value. Based on all of the above information, it can 
be seen that Re approaches the value of R, so tip losses are essentially neglected. 
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The third consideration is ground effect. It is known that the induced velocity and thrust required 
to hover are considerably reduced in ground effect. The ratio of induced velocity to that which would have 
occurred in free air is shown in figure 1 1 as a function of the radial position and the ratio of rotor height to 
rotor radius. The ratio of thrust in ground effect to thrust in free air at a given free air power setting as a 
function of rotor height and thrust coefficient/solidity is also shown in figure 1 1 (ref 3) The deter- 
mination of induced velocity ratio is subject to the assumption that induced velocity is considered to be 
uniform for the purposes of this study. Induced velocity, vj, is actually directly proportional to radial 
ocation (e.g., yi - kr, from r - 0 to r - R) such that the average value is located at r/R = 2/3. This 
average value is used as the uniform value, and the induced velocity ratio can be obtained from figure 1 1 
based on r/R = 2/3. The rotor height, h, at this location is 22 ft based on the fact that the hovering height 
of the lowest point of the Da Vinci II, h 0 , is 10 ft off the ground. Therefore, h/r = 0.5, and the induced 
velocity ratio vAi .obtained from figure 1 1 is approximately 0.35. The deteimination of thrust ratio is a 
function of h/r and CjJc. The value of C T Ja is 0. 17 and the thrust ratio at a given free-air power 
setting obtained from figure 1 1 is approximately 1.2. The inverse of this value represents the ratio of 
trust required to hover in ground effect to thrust required to hover in free air, and is approximately 0.833. 
This constitutes a 16.7% decrease in the thrust required and the thrust coefficient when in ground effect. 

i fourth consid eration is the effect of the control surfaces on thrust and induced velocity The 
results of this report show that each control surface area should be 1 1.5 ft 2 in order to obtain a stable well 
behaved system The control surfaces are set at the zero-lift angle of attack when the rotors are leveled or 
when the control system is off. The control surfaces move differentially when actuated, therefore the total 
lift generated equals zero at all times. Drag is generated and must be included in the profile drag calcula- 
tion used above for the development of thrust required to hover and induced velocity. Equation (25) cal- 
culates this value of drag, as follows: 4 K } 


D = D 0 + D cs (max) 

where D cs (max) = c,j 0 p[D(R + l cs /2)] 2 S cs 

S cs = 1 1.5 ft 2 
lcs = Sqs/C = 3.83 ft 
D cs (max) = 0.51 lb 


This value comprises roughly 6.0% of D 0 . Based on this information 
surfaces to thrust and induced velocity are neglected. 


the contributions of the control 


, /^fore, Ihe expressions and values for l x , I y , T, Or. and v* based on aeroelasnci.y and ground 
effect for hovenng flight are restated by equations (26) to (30), as follows: * S 


- 221.0 + 1 133.3(sin \j/) 2 + 7253.2(cos y) 2 

(26) 

= 221.0 + 1 133.3(cos y) 2 + 7253.2(sin y)2 

(27) 

= 290.47/1.2 = 242.06 lb 

(28) 

= 0.00489/1.2 = 0.00408 

(29) 

= 0.35 x vj^ = 0.69 ft/sec 

(30) 

where v ioo = 2Vf 0< /(pA)/3 = 1.96 
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APPENDIX C 


AIRCRAFT PARAMETER AND DERIVATIVE VALUES 


This appendix contains two tables. The first table, table 3, lists values for all aircraft specifications 
of the Da Vinci II. The second table, table 4, lists values for all stability derivatives, as well as control 
system parameters used during the course of this study. 


TABLE 3.- AIRCRAFT SPECIFICATIONS 


Symbol 

Value 

Symbol 

Value 

a 

6.45 

h 

22.0 

A 

14,102.19 

ho 

10.0 

hR 

1.5 

m 

8.85 

ib 

2323.5 

mb 

5.12 

p 

0.17452 

nv 

1.55 

c 

3.0 

m t 

0.31 

Cdo 

0.01 

Q 

0.6283 

Ct 

0.00408 

N 

2 

CToo 

0.00489 

a 

0.0285 

C.F. 

28.75 

R 

67.0 

d 

5.5 

P 

0.002378 

d a 

4.0 


-0.0164 

db 

1.5 


0.1745 

dc 

12.0 

W 

285.00 

dp 

6.0 

tmr 

0.33 

Do 

8.5 

T 

242.06 

e 

1.0 

Too 

290.47 

y 

399.07 

Vi 

0.69 

g 

32.2 

Vioo 

1.96 
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TABLE 4.- STABILITY DERIVATIVES AND 
CONTROL SYSTEM PARAMETERS 


Symbol 

Value 

^cs 

les 

Del x /sin \| / 

D^Iy/cos y 

Ka 

K b 

Kc,iIx/(cos \j/) 2 
K c , 2 lx/(sin \) / cos y) 
Kd,lly/(sin \j/ cos \j/) 
Kd,2ly/(sin \j /) 2 
Kpl x 
L v Ix 

Mqly 

M u r y 

N r 

N v 

S 

vt 

Zw 

11.5 

3.83 

22790.0 
-22790.0 

0.017452 

-0.017452 

-53296.0 

53296.0 

53296.0 

-53296.0 

-149.0 

-988.6 

-149.0 

988.6 

0 

0 

1.70 

-11.81 

-1.70 

-11.81 

-2.92 
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APPENDIX D 


CHARACTERISTIC POLYNOMIAL DEVELOPMENT 


The development of the characteristic polynomials for the longitudinal and lateral modes of the 
unaugmented and augmented configurations of the Da Vinci II are described in this appendix. The equa- 
tions describing the kinematics and control system were rearranged in state form such that the characteristic 
polynomials are defined by {det fsl - A]} =0 (ref. 7). 


The longitudinal mode of the unaugmented configuration is considered first, as described by equa- 
tions (31) to (34). 


ii = X u u - g0 + X q q + X u u' 


(31) 


q = M u u + M q q + D 08 cs,d + M u u 


9 =q 


xi 

X2 

_x 3 _ 


X u 

0 

M u 



0 

0 

De 



u 


" 1 0 o ” 


xi 

e 

= 

0 1 0 


X2 

_e_ 


_ 0 0 1 _ 


_X3_ 


Solving for {det[sl - A]) = 0: 


{det[sl - A] } = s(s - X u )(s - M q ) - sM u X q + gM u = 0 


(32) 

(33) 


= 1 - M u 


(* q s - g) 

s(s - X u )(s - M q ) 


= 0 


(34) 


The lateral mode characteristic polynomial of the unaugmented configuration was developed in a 
similar fashion and is described by equation (35): 


{ det[ si - A] } = s(s - Y v )(s - Lp) - sL v Y p + gL v — 0 


_i T ( y n s + g> o 

- 1 " Lv s(s- Y v )(s- L p ) 


(35) 


The augmented configuration is considered next. The equations describing the kinematics and 
control system are combined and the longitudinal and lateral modes become coupled, as described by 
equations (36) to (40), as follows: 
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XI 
X2 
X3 
X 4 
X5 
ux 6 J 

XI 
x 2 
x 3 
X4 
X5 
1— X6_J 


U “ x u u - g0 + Xq0 + X u u’ 

0 = M u u + Mq0 + Kdj*]) + K d t2 0 + D05 CSjd + M u u' 
v = Y v v + g<{» + Y p <|) 

<j) = L v v + Lp0 + K c> i<{) + Kc )2 0 + D<j)8 CSj{ i 

X u 0 
0 Y v 
0 0 


-g 

0 

0 


Xq 

0 

1 


0 0 


g 

0 


0 


M u 0 Kd >2 Mq Kd,i 0 
0 0 0 0 0 1 

L 0 Ly K c 2 0 K c i Lp _] 


XI 


” X u 

0 

x 2 


0 

0 

x 3 


0 

0 

X4 


M u 

E>0 

x 5 


0 

0 

-X6_ 


_ 0 

D<{> _ 


Locs.dJ 


1 0 0 0 0 0 

0 1 0 0 0 0 

0 0 1 0 0 0 

0 0 0 1 0 0 

0 0 0 0 1 0 

L 0 0 0 0 0 1 JL<j)J 


u 

V 

0 

0 

<J> 


(36) 

(37) 

(38) 

(39) 


Solving for {det[sl - A] } = 0 using pivotal condensation (also called the method of Chio) for the 
augmented configuration yields the following characteristic polynomial (ref. 4): 

1 -— K c , 2 K d ,i(s-X„)(s- Y v ) 

((s - X u )[s(s - M q ) - K d ,2l - M„(X q s - g)} { (s - Y v )[s(s - L p ) - K Ct i] - L v (Y p s + gj] “ 0 


(40) 
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APPENDIX E 


DISCRETE SIMULATION DESCRIPTION 


A discrete simulation was developed at the NASA Ames Research Center Flight Systems and 
Simulation Research Division to study the stability characteristics of the augmented configuration of the 
Da Vinci II and to size control system design parameters in order to achieve a stable, well behaved 
system. 

The simulation was written in Fortran and was comprised of one main program and four subrou- 
tines called by the main program. The main program called the four subroutines in the same order in 
12 separate loops. The function of each subroutine and the purpose of the loops are described in this 
appendix; a flow diagram is shown in figure 12. 

The function of the first subroutine is to introduce perturbations to the system of equations 
describing the kinematic model. The perturbations are in the form of body axis velocity wind gusts, u' 
and w'. \|/ is integrated in this subroutine so that at a particular point in time of a simulation run, a value 
of u' or w' is introduced. This value is reset back to zero in a finite amount of time so that the input 
resembles a step. 

The function of the second subroutine is to model the control system of the Da Vinci II. This sub- 
routine calculates rotor tip height, rotor tip height difference, control surface deflections, and restoring- 
moments and associated accelerations (pfl> and qfb). The control surface deflection that is calculated from 
the rotor height difference is differentiated so that a rate limit can be imposed on the system. The rate- 
limited value is then integrated so that a position limit can be imposed. These two parameters, along with 
control surface area, were varied in order to define design specifications of the actuators necessary to 
achieve the desired, stable behavior. 

The function of the third subroutine is to model the kinematic equations. This subroutine calculates 
body axis accelerations and uses Adam's-Bashforth integration to obtain body axis rates and positions 
(ref. 6). Some of these rates and positions are used by the kinematic model and control system for each 
successive cycle. The cycle time used for this simulation is 0.05 sec. 

The function of the fourth subroutine is to print the aircraft states and stability derivatives of the Da 
Vinci II. This subroutine prints a line each second stating that the control system is on, if that be the case. 
Finally, this subroutine will print a "CRASH" message if a rotor hits the ground, and will stop printing 
after that printout. However, the simulation continues to cycle until the run is complete. 

These subroutines are called in 12 separate loops. The first loop initializes arrays and variables in 
each subroutine as well as the main program. The second loop sets appropriate variables to their trim val- 
ues for hovering flight. The trim conditions for this simulation are as follows: 

h 0 = 10.0 ft 

Q = 0.6283 rad/sec 

Body axis accelerations, rates, and positions are initialized to zero 
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Variables used in numerical integration and differentiation are initialized to zero 
The control system is turned off 

The third loop introduces the u' perturbation to the kinematic model and allows the simulation to cycle for 
1 min, printing aircraft states and stability derivatives each second. These three loops are repeated using 
the w' perturbation. The whole process is then repeated with the control system turned on. 

The actual Fortran code for this simulation is contained herein. Fortran variables, definitions, and 
associated units are defined in the Comments sections of the main program and four subroutines. Com- 
ments are also included in the code in order to explain the logic. 
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uuuouuoouuuuuuuouuuuuouuuououuuuuouuuuuuuuuuouu 


TITLE 

PROGRAM COMMAND 


CREATION/MODIFICATION LOG: 
DATE NAME REMARKS 

3/88 J. TOTAH (NASA) WRITTEN 


INTRODUCTION: 


THIS PROGRAM CONTROLS THE CALLING SEQUENCE OF THE SUBROUTINES 
THAT COMPRISE THE MODEL FOR THE CALIFORNIA POLYTECHNIC HUMAN- 
POWERED HELICOPTER PROJECT, THE DAVINCI II. THE CALLING 
SEQUENCE IS AS FOLLOWS: 

PROGRAM SUBROUTINES 


COMMAND - CHKDYN 

- CONTR2 

- AER02 
- PRTOUT 

THE SUBROUTINES' FUNCTIONS ARE DEFINED AS FOLLOWS: 

CHKDYN: CALCULATES CONTROL SURFACE PERTURBATIONS 
CONTR2 : CALCULATES RESTORING MOMENT AND ACCELERATION FEEDBACK 
AER02 : CALCULATES AIRCRAFT STATES, BODY AXIS AND EARTH AXIS 
PRTOUT : PRINTS AIRCRAFT STATES AND STABILITY DERIVATIVES 

DEFINITIONS OF VARIABLES ARE AS FOLLOWS: 


INPUTS: 

ITRIM I A (06) NUMBER OF CYCLES TO TRIM CYC DATA 

I NTT IA(07) NUMBER OF CYCLES TO INITIALIZE CYC DATA 

ID IA(05) SWITCH TO STOP PRTOUT UPON IMPACT N/A PRTOUT 


RPM DA(033) ROTOR SPEED 

DT2 DA(088) CYCLE TIME 


RAD/SEC CONTR2 
SEC CONTR2 


OUTPUTS: 


CIMODE IA(01) 

C IPART IA(02) 

C ICHKDYN IA(03) 

C ISAS IA(08) 

CIDEL1 IA(13) 


MODE CTRL. INT. (-:IC 0:HLD +:OP) 
COUNTER USED IN PRTOUT 
SWITCH TO ACTIVATE 'CHKDYN' 
SWITCH TO TURN SAS ON AND OFF 
LONGITUDINAL WIND GUST SWITCH 


N/A COMMAND 
N/A COMMAND 
N/A COMMAND 
N/A COMMAND 
N/A COMMAND 
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CIDEL2 IA( 14) VERTICAL WIND GUST SWITCH N/A 

C 

C LOCALS: 

C 

C IDYNAMIC IA(04) NUMBER OF CYCLES FOR 6 REVOLUTION CYC 
C 
C 
C 

COMMON /DA VINCI/DA( 1 50) 

COMMON /IFIXED/IA(20) 

C 

c 

c 

EQUIVALENCE (IA(Ol).IMODE ) 

EQUIVALENCE (IA(02),IPART ) 

EQUIVALENCE (IA(03),ICHKDYN ) 

EQUIVALENCE (IA(05),ID ) 

EQUIVALENCE (IA(08),ISAS ) 

EQUIVALENCE (IA(13),IDEL1 ) 

EQUIVALENCE (IA(14),IDEL2 ) 

C 

EQUIVALENCE (DA(033),RPM ) 

EQUIVALENCE (DA(088),DT2 ) 

C 

C 

DATA INIT, ITRIM / 10 , 10/ 

DATA DA ,IA /1 50*0., 10*0./ 

C 

C**** BEGIN EXECUTABLE CODE 
C 

IDYNAMIC = ((12.*3.1415)/RPM) * (1./DT2) 

C 

C 

ISAS = 0 

DO 25 K= 1,2 

IF (K .EQ. 2) ISAS = 1 

IDEL1 = 1 

IDEL2 = 0 

DO 20 J= 1,2 

IF (J .EQ. 2) IDEL1 = 0; IDEL2 = 1 
C 

c**** CYCLE IN I.C. MODE TO INITIALIZE FILTERS AND VARIABLES 
C 

IMODE =-l 

ICHKDYN = 0 

IPART = 0 

C 

DO 5 I = 1,INIT 

IPART = IPART + 1 

CALL CHKDYN(I) 

CALL CONTR2 


COMMAND 

COMMAND 
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ooooo onn 


CALL AER02 
CALL PRTOUT 
CONTINUE 


c**** trim the aircraft 

L IMODE = -1 

ICHKDYN = 0 
IP ART = 0 


10 

C 

£**** 


DO 10 1 = ™ 

TP ART = IPART 

CALL CHKDYN(I) 
CALL CONTR2 
CALL AER02 
CALL PRTOUT 
CONTINUE 


+ 1 


perform dynamic checks on the aircraft 


IMODE = ; 

ICHKDYN = 1 

IPART = -1 

ID =0 


c 


14 

15 

20 

25 


DO 15 1 

IF(ID.EQT)G0T0 14 

TP ART = IPART + 1 
CALL CHKDYN(I) 

CALL CONTR2 
CALL AER02 
CALL PRTOUT 
CONTINUE 
CONTINUE 


CONTINUE 

CONTINUE 


STOP 
END 

TITLE 

SUBROUTINE CHKDYN 

SUBROUTINE CHKDYN(I) 

c' creation/modification LOG 
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c DATE name remarks 

C— - 3 - 8 - 8 - J ' T0TAH (NASA) WRITTEN 

c - — 

c INTRODUCTION: 

C W ™ LON OmjDINAL AND 

c d ™ons of variables are as follows- 

c INPUTS: 

ClD^i D yA(13) (04) LONGmrr? IN ™ 0D UCE PERTURBATIONS N/A C0MMa ND 

?““■*<“> asasisa™ ssa 

CGT7ME DA(1I9) beginning of input WA command 

c OUTPUTS: SEC DAT A 

Crem Slif! <!r?r n - CHKDYN 

?s S Sssasli B ass 

c DA018) vertical wind GUST US CHKDYN 

c LOCALS: 1/2 CHKDYN 

c NTIME IA(09) TlA/rn 

CNNIME IA(10) ^AT,ONT, E N?r NGOFINP ^ CVC CHKDVN 
C CYC CHKDYN 

c sssar* 

c 

c 


c 


I&e E ^»S DYN J 

EQUIVALENCE (IA(H),IDEL2 j 


EQUIVALENCE (DAfO'^R'i r 
EQUIVALENCE (DA 032 pi 
EQUIVALENCE (DA(03 R i 
EQUIV ALENCE (DA(038) X 
EQUIVALENCE DA (039 X 
equivalence (DA(040 V 
equivalence DA(04? ’? 
EQUIVALENT ns nfi ’B 



) 

) 

) 

) 

) 

) 

) 

) 

) 
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EQUIVALENCE (DA(049),XXI ) 

EQUIVALENCE (DA(050),YYI ) 

EQUIVALENCE (DA(088),DT2 ) 

EQUIVALENCE (DA(089),PSIR ) 

EQUIVALENCE (DA(095),CSAREA ) 

EQUIVALENCE (DA(096),XLCS ) 

EQUIVALENCE (DA(097),CG ) 

EQUIVALENCE (DA(1 1 1),AD ) 

EQUIVALENCE (DA(112),BD ) 

EQUIVALENCE (DA(1 1 3), CD ) 

EQUIVALENCE (DA(1 17), DELI ) 

EQUIVALENCE (DA(1 1 8),DEL2 ) 


DATAGTIME/ 1.25/ 

C**** jnertIA CALCULATIONS AS A FUNCTION OF PSI 
C 

PSIR = PSIR + RPM*DT2 

1 ^^XMR^^^DIUS*COS^BETA)* COS(PSnl))* *2)/3 . + (AD+CD/2.-CG)**2 
1 2. amk l RA DIUS*COS(BETA)*SIN(PSIR))**2)/12.) + 

3 2.*XMT*(((RADIUS*COS(BETA)*COS(PSIR))**2) + (AD+CD-CG)**2) 

1 ^^*XN^^WRADIU^:CO^(^ETAy |e SD^PSIR))**2)/3. + (AD+CD/2.-CG)**2 

2 4 ^ RADIUS*COS(BETA)*COS(PSIR))**2)/12^) + 

3 2.*XMT*(((RADIUS*COS(BETA)*SIN(PSIR))**2) + (AD+CD-CG)**2) 

c 

IF (ICHKDYN .NE. 1) GO TO 10 

C**** §eT TIME OF INPUT AND AMPLITUDE OF STEP PURTURBATIONS 
C 

NTIME = GTIME/DT2 
NNIME = NTIME* 3. 

C 

C **=K* DYNAMIC check 

c 

IF (I .LE. NTIME) GO TO 10 
IF (I .GE. NNIME) DELI = 0.; DEL2 = 0.; GO TO 10 
r 

c **** XHE VALUE 7.333 FT/S REFERS TO A 5 MPH WIND 
C 

DELI = IDEL1 *7.333 
DEL2 = IDEL2*7.333 
C 

10 CONTINUE 
C 

c 

c 


RETURN 
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END 


C 

C 

C 

c 

c 

c 

c 


TITLE 

SUBROUTINE CONTR2 


SUBROUTINE CONTR2 



C CREATION/MODIFICATION LOG: 

C DATE NAME REMARKS 

C 3/88 J. TOTAH WRITTEN 

C 

C INTRODUCTION: 

C 



C DEFINITIONS OF VARIABLES ARE AS FOLLOWS: 

C INPUTS: 

C 


CISAS 

C 

C RADIUS 
CN 

C CHORD 

CCDO 

CRHO 

CRPM 

C WAIT 

C WROTOR 

CWBODY 

C WTIP 

C BETA 

CGEF1 

CGEF2 

CTHETO 

C DCLDA 

CE 

CXXI 

C YYI 

CG 

CTHETR 
CPHIR 
CDT2 
CPSIR 
C CSAREA 


IA(07) AUGMENTATION ON/OFF SWITCH 


DA(028) 

DA(029) 

DA(030) 

DA(031) 

DA(032) 

DA(033) 

DA(034) 

DA(035) 

DA(036) 

DA(037) 

DA(042) 

DA(043) 

DA(044) 

DA(045) 

DA(046) 

DA(047) 

DA(049) 

DA(050) 

DA(051) 

DA(072) 

DA(073) 

DA(088) 

DA(089) 

DA(095) 


MAIN ROTOR RADIUS 
NUMBER OF BLADES 
MAIN ROTOR CHORD 


PKUHLE DRAG COEFFICIENT 
AIR DENSITY AT SEA LEVEL 
ROTOR SPEED 


HELICOPTER TOTAL WEIGHT AT HOVE 
ROTOR BLADE WEIGHT 
PILOT AND FRAME WEIGHT 
WEIGHT AT ROTOR TIP 
CONING ANGLE 
VI GROUND EFFECT FACTOR 
THRUST GROUND EFFECT FACTOR 
MAIN ROTOR INITIAL PITCH ANGLE 
ROTOR LIFT-CUR VE-SLOPE 
HINGE OFFSET 


ROLL MASS MOMENT-OF-INERTIA 
PITCH MASS MOMENT-OF-INERTIA 
GRAVITY 


PITCH ANGLE, BODY AXIS 
ROLL ANGLE, BODY AXIS 
CYCLE TIME 


ROTOR POSITION IN THE TPP 
CONTROL SURFACE AREA 


N/A 

DATA 

FT 

DATA 

N/A 

DATA 

FT 

DATA 

N/A 

DATA 

SLUGS/FT3 

DATA 

FT/SEC 

DATA 

LB 

DATA 

LB 

DATA 

LB 

DATA 

LB 

DATA 

RAD 

DATA 

N/A 

DATA 

N/A 

DATA 

RAD 

DATA 

1/RAD 

DATA 

N/A 

DATA 

FT2-SLUGS 

CHKDYN 

FT2-SLUGS 

CHKDYN 

FT/S 2 

DATA 

RAD 

AER02 

RAD 

AER02 

SEC 

DATA 

RAD 

CHKDYN 

FT2 

DATA 
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c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


GRL DA(108) ACTUATOR RATE LIMIT 

GPL DA(109) ACTUATOR POSITION LIMIT 

GKA DA(110) HEIGHT-TO- ANGLE GAIN 

AD DA(lll) DISTANCE FROM BASE TO HUB 

BD DA(1 12) DISTANCE FROM BASE TO PILOT CG 

CD DA(1 13) DISTANCE FROM HUB TO ROTOR TIP 

PSIRIC DA(1 14) INITIAL ROTOR POSITION IN TPP 

DELI DA(1 17) LONGITUDINAL WIND GUST 

DEL2 DA(118) VERTICAL WIND GUST 


OUTPUTS: 


PBDFB 

QBDFB 


DA(093) ROLL ACCELERATION FEEDBACK 
DA(094) PITCH ACCELERATION FEEDBACK 


LOCALS: 

XLCS DA(096) LENGTH OF CONTROL SURFACE 

CG DA(097) HELICOPTER CG 

HITEA DA(098) SIGNAL TO ACTUATOR 

DHITE DA(099) DIFFERENTIATED ACTUATOR SIGNAL 

HITE DA 100) CONTROL SURFACE POSITION 

RREST DA(104) RESTORING FORCE OF CONTROLS 

PM DA115 ROLL MOMENT DUE TO CONTROLS 

QM DA(1 16) PITCH MOMENT DUE TO CONTROLS 


COMMON /IFIXED/IA(20) 
COMMON /D A V IN CI/D A ( 1 50) 


RAD/S 

RAD 

RAD/FT 

FT 

FT 

FT 

RAD 

FT/S 

FT/S 


RAD/S 2 
RAD/S2 


FT 

FT 

RAD 

RAD/S 

RAD 

LB 

FT-LB 

FT-LB 


DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

DATA 

CHKDYN 

CHKDYN 


CONTR2 

CONTR2 


CONTR2 

CONTR2 

CONTR2 

CONTR2 

CONTR2 

CONTR2 

CONTR2 

CONTR2 


EQUIVALENCE (IA(Ol)JMODE ) 

EQUIVALENCE (IA(08),ISAS ) 

EQUIVALENCE (IA(13),IDEL1 ) 

EQUIVALENCE (IA(14),IDEL2 ) 

EQUIVALENCE (DA(028), RADIUS ) 

EQUIVALENCE (DA(029),N ) 

EQUIVALENCE (DA(030), CHORD ) 

EQUIVALENCE (DA(031),CDO ) 

EQUIVALENCE (DA(032),RHO ) 

EQUIVALENCE (DA(033),RPM ) 

EQUIVALENCE (DA(034)WAIT ) 

EQUIVALENCE (DA(038),XMASS ) 

EQUIVALENCE (DA(039),XMR ) 

EQUIVALENCE (DA(040),XMB ) 

EQUIVALENCE (DA(041),XMT ) 

EQUIVALENCE (DA(042),BETA ) 

EQUIVALENCE (DA(043),GEF1 ) 

EQUIVALENCE (DA(044),GEF2 ) 
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ODD 


C 

c 

c 


EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 

EQUIVALENCE 


(DA(045),THETO 

(DA(046),DCLDA 

(DA(047),E 

(DA(048),HR 

(DA(049),XXI 

(DA(050),YYI 

(DA(051),G 

(DA (07 2) ,THETR 

(DA(073),PHIR 

(DA(088),DT2 

(DA(089),PSIR 

(DA(093),PBDFB 

(DA(094),QBDFB 

(DA(095),CSAREA 

(DA(096),XLCS 

(DA(097),CG 

(DA( 100), HITE 

(DA(111),AD 

(DA(112),BD 

(DA(1 13), CD 

(DA(1 17), DELI 

(DA(118),DEL2 


) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 


DATA CDO 
DATAE 
DATA THETO 
DATA PSIRIC 
DATA WBODY 
DATA WTTP 
DATA DCLDA 
DATA CSAREA 
DATA DT2 , 
DATAG 
DATA BD 
DATA GPL 
DATA ISAS 
DATA BETA 


, N 

/ 0.01 

, GEF1 

/ 1.0 

, GEF2 

/ 0.1745 

, GKA 

/ -0.031415 

, WROTOR 

/1 65.0 

, WAIT 

/ 10.0 

, RADIUS 

/ 6.45 

, RPM 

/ 11.5 

RHO 

/ 0.05 

, AD 

/ 32.2 

, CD 

/ 1.5 

, GRL 

/ 0.04 

, CHORD 

/ 0 

/ 0.1745/ 


. 2 / 

, 0.35 / 

, 1.2 / 

, -0.01745 / 
, 50.0 / 

,285.0 / 

, 67.0 / 

, 0.6283 / 

, 0.002378 / 

, 4.0 / 

, 12.0 / 

, 0.24 / 

, 3.0 / 


IC MODE — 


IF (IMODE) 100,300,200 

100 CONTINUE 
C 

C**** FILTER INITIALIZATION 

p SIR = PSIRIC 

HITEA = 0 
HITEAP = 0 
C 
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HITE = 0. 

DHITE = 0. 

DHITEP = DHITE 
C 

c **** DESIGN specification calculations 

= SQRT(CSAREA) 

= WBODY / G 
= WROTOR/ G 
= WTIP /G 
_ v/AIT / G 

= 2.*(XMR*(AD+CD/2.) + XMT*(AD+CD))/XMASS + 
XMB*BD/XMASS 
= ABS(CG-AD) 

— OPERATE MODE 


XLCS 

XMB 

XMR 

XMT 

XMASS 

CG 

1 

HR 


200 CONTINUE 

£.«, CALCULATION OF DIFFERENTIAL ROTOR HEIGHT COMMAND TO ACTUATORS 


c 


HITEA 
1 *ISAS 


2.*RADIUS*(PHIR*COS(PSIR) - THETR*SIN(PSIR))*GKA 


C**** rate AND POSITION LIMIT CALCULATIONS 
C DHITE = (HITEA - HITEAP)/DT2 

H^ABSIDHITCLOE. GRL) DHITE = SIGN(GRL,DHITE) 


c 


HITE 

DHITEP 


= HITE + 0.5*DT2*(3.*DHITE - DHITEP) 
= DHITE 


IF (ABS(HITE) .GE. GPL) HITE = SIGN(GPL,HITE) 


£**** CALC. OF FEEDBACK ACCELERATIONS FROM CONTROL SURFACE DEFLECTIONS 

C RREST = DCLDA*HITE*RHO*((RPM*(RADIUS+XLCS/2.))**2.)*CSAREA 

C OM -RREST*(RADIUS+XLCS/2.)*(SIN(PSIR)) 

PM = RREST* (RADIUS+XLCS/2.)* (COS (PSIR)) 


QBDFB = IDEL1 *QM/YYI 

PBDFB = IDEL1*PM/XXI 


C 

C 

C 

C 


300 


CONTINUE 

RETURN 

END 


C 

C 
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c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

C- 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 


TITLE 

SUBROUTINE AER02 


SUBROUTINE AER02 


CREATION/MODIFICATION LOG; 
DATE NAME REMARKS 
3/88 J. TOTAH WRITTEN 


INTRODUCTION: 

iSKS FORCES, AND 

DEFINITIONS OF VARIABLES ARE AS FOLLOWS: 

INPUTS: 


PBDFB 

QBDFB 


OUTPUTS: 


DA(093) ROLL ACCELERATION FEEDBACK 
DA(094) PITCH ACCELERATION HEEDBA(^ 


SDXU 

SDZW 

SDMU 

SDMQ 

SDXQ 

SDYV 

SDYP 

SDLV 

SDLP 

UBD 

VBD 

WBD 

PBD 

QBD 

UB 

VB 

WB 

PB 

QB 

THETR 

PHIR 

XE 

YE 


DA (001) 

DA(002) 

DA(003) 

DA (004) 

DA(005) 

DA(006) 

DA(007) 

DA (008) 

DA (009) 

DA(052) 

DA(053) 

DA(054) 

DA(055) 

DA(056) 

DA (057) 

DA(058) 

DA(059) 

DA(060) 

DA(061) 

DA(072) 

DA(073) 

DA(080) 

DA(081) 


LONG. DRAG DAMPING 
VERTICAL DAMPING 

KSo VELOC ^ STAB[LrrY 

PARTIAL X-FORCE WRT PITCH RATE 
LATERAL DRAG DAMPING 
PARTIAL Y-FORCE WRT ROLL RATE 
LATERAL VELOCITY STABILITY 
ROLL DAMPING 
FORWARD ACCELERATION 
SIDE ACCELERATION 
VERTICAL ACCELERATION 
ROLL ACCELERATION 
PITCH ACCELERATION 
FWD BODY VELOCITY 
LFT BODY VELOCITY 
DWN BODY VELOCITY 
ROLL RATE 
PITCH RATE 
PITCH ANGLE 
ROLL ANGLE 
LONGITUDINAL POSITION 
LATERAL POSITION 


BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

BODY 

EARTH 

EARTH 



CONTR2 


CONTR2 

RAD/SEC 

AER02 

RAD/SEC 

AER02 

RAD/SEC/FT 

AER02 

1/RAD/SEC 

AER02 

ft/sec/rad 

AER02 

RAD/SEC 

AER02 

FT/SEC/RAD 

AER02 

RAD/S/FT 

AER02 

1/RAD/SEC 

AER02 

AXIS R/S2 

AER02 

AXIS FT/S2 

aero: 

AXIS FT/S2 

AER02 

AXIS R/S2 

AER02 

AXIS R/S2 

AER02 

AXIS FT/S 

AER02 

AXIS FT/S 

AER02 

AXIS FT/S 

AER02 

AXIS RAD/S 

AER02 

AXIS RAD/S 

AER02 

AXIS RAD 

AER02 

AXIS RAD 

AERD2 

AXIS FT 

AER02 

AXIS Fi’ 

AER02 
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CZE 
C FIX 
CFTY 
CFTZ 
C TIL 
CTTM 
C 

C LOCALS 
C 

CDADU 
CDHDU 
CDADQ 
CDHDQ 
CDCTDW 
CDMDU 
CDMDQ 
CBBI 
C GLOCH 
CCT 
CCF 
CXED 
CYED 
CZED 
C 


DA(082) 
DA(083) 
DA(084) 
DA(085) 
DA(086) 
DA (087) 


VERTICAL POSITION 

AERODYNAMIC X-FORCE 

AERODYNAMIC SIDE FORCE 
AERODYNAMIC VERTICAL FORCE 
AERODYNAMIC ROLL MOMENT 
AERODYNAMIC PITCH MOMENT 


EARTH AXIS FT 
BODY AXIS LBF 
BODY AXIS LBF 
BODY AXIS LBF 
BODY AXIS FT-LBF 
BODY AXIS FT-LBF 


C 

C 


DA(010) 

DA(011) 

DA(012) 

DA(013) 

DA(014) 

DA(015) 

DA(016) 

DA(024) 

DA(025) 

DA(026) 

DA (027) 

DA(074) 

DA(075) 

DA(076) 


« OE FRFORCE DUE TO UB **£*** 

FLAPPING COEFF. FORCE DUE TO QB u( ||^s^ ad 
HUB FORCE DUE TO QB SLUGS-FT/S/RAU 

THRUST COEFFICIENT DUE TO WB 
SHAFT MOMENT DUE TO UB 

ROTOR BLADE INERTIA 
LOCH NUMBER 
COEFFICIENT OF THRUST 
CENTRIFUGAL FORCE 
LONGITUDINAL SPEED 

lateral speed 

VERTICAL SPEED 


SLUGS-FT-RAD/S 

N/A 

FT2-SLUGS 

N/A 

N/A 

LB 

AXIS FT/S 
AXIS FT/S 
AXIS FT/S 


EARTH 

EARTH 

EARTH 


COMMON /IFIXED/I A (20) 
COMMON /DAVINCI/DA(150) 


C 

C 

C 

C 


EQUIVALENCE (IA(OOl)JMODE 



) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 

) 


AER02 

AER02 

AER02 

AER02 

AER02 

AER02 


AER02 

AER02 

AER02 

AER02 

AER02 

AER02 

AER02 

AER02 

AER02 

AER02 

AER02 

AER02 

AER02 

AER02 
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(DA(045),THETO 

(DA(046),DCLDA 

(DA(047),E 

(DA(048),HR 

(DA (049), XXI 

(DA(050),YYI 


EQUIVALENCE 

equivalence 
equivalence 
equivalence 

EQUIVALENCE v-^vuju) r 
EQUIVALENCE (DA(051)G 

pffir L ^ NCE (DAIOsiiuBD 

EQUIVALENCE (DA(053) VBD 
EQUIVALENCE - • V J ’ 

equivalence 
equivalence 
equivalence 

EQUIVALENCE vh 

EQUIVALENCE (DA(059) WB 

EQUIVALENCE (DAWoJiPB 

EQUIVALENCE (DA(061)OB 

(5 A(072 *'thetr 


(DA(054),WBD 

(DA(055),PBD 

(DA(056),QBD 

(DA(057),UB 

(DA(058),VB 


EQUIVALENCE (DA(073)PHIR 

f^,^ NCE ( da S™ 

EQUIVALENCE (DA(081)YE 
EQUIVALENCE (DA 082) ZE 
EQUIVALENCE (DA(083) FTX 
EQUIVALENCE (DA (084) FT Y 

S3KSSS SSSIS 
isiKtmg 

equivalence 
equivalence 


(DA(088),DT2 
(DA(089),PSIR 
(DA(093),PBDFB 


EQUIVALENCE (DA 094) QBDFB 
EQUIVALENCE (DA(1 nXDELl 
EQUIVALENCE (DA(1 18),DEL2 
EQUIVALENCE (DA(123 )|ho 


IC MODE 


c IF (IMODE) 100,300,200 

100 CONTINUE 

C**** FILTER INITIALIZATIONS 

UBD = 0 
VBD = o' 

WBD = 0 
PBD = o' 

PB = 0. 

Qbd = o. 

QB =0. 
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c 

UBDP = UBD 
VBDP = VBD 
WBDP = WBD 
PBDP = PBD 
PBP = PB 
QBDP = QBD 
QBP = QB 
C 

UB = 0. 

WB = 0. 

VB = 0. 

PB = 0. 

PHIR = 0. 

QB = 0. 

THETR =0. 

C 

XED = 0. 

YED = 0. 

ZED = 0. 

C 

XEDP = XED 
YEDP = YED 
ZEDP = ZED 
C 

XE = 0. 

YE =0. 

ZE = HO 
C 

c**** STABILITY DERIVATIVE CALCULATIONS 
C 

AREA = 3.1415*(RADIUS**2) 

SIGMA = N*CHORD/(3. 1415*RADIUS) 

DO = CDO*RHO*CHORD*RADIUS*((RPM*RADIUS)**2) 

TOO = S QRT ( W AIT* * 2 + DO**2)/COS(BETA) 

THRUST = TOO/GEF2 

VI = GEFl*2.*SQRT(TOO/(RHO*AREA))/3. 

XIN = -VI/(RPM*RADIUS) 

DADU = 8.*THETO/3. + 2.*XIN 

DHDU = RHO*SIGMA*AREA*RPM*RADIUS*CDO/4. 

SDXU = -(THRUST*DADU + DHDU)/XMASS 


BBI = XMR*(RADIUS**2)/3. 

GLOCH = RHO*DCLDA*CHORD*(RADIUS**4)/BBI 

DADQ = -16./(GLOCH*RPM) 

DHDQ = -RHO*DCLDA*AREA*SIGMA*((RPM*RADIUS)**2)*XIN/ 

1 (2.*GLOCH*RPM) 
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SDXQ = -(THRUST*DADQ + DHDQ)/XMASS 
SDYV = SDXU 

SDYP = -SDXQ 

CT = TOO/(RHO*AREA*((RPM*RADIUS)**2))/GEF2 

DCTDW = DCLDA*SIGMA/(8. + DCLDA*SIGMA*SQRT(2./CT)/2.) 

SDZW = -RHO*AREA*RPM*RADIUS*DCTDW/XMASS 

CF = (XMR/2. + XMT)*RADIUS*(RPM**2) 

DMDU = N*E*RADIUS*CF*DADU/2. 

SDMU = (HR*(DHDU + THRUST*DADU) + DMDU)/YYI 

DMDQ = -8.*N*E*RADIUS*CF/(GLOCH*RPM) 

SDMQ = (HR*(DHDQ + THRUST*DADQ) + DMDQ)/YYI 

SDLV = -SDMU*YYI/XXI 

SDLP = SDMQ*YY1/XXI 

OPERATE MODE 

200 CONTINUE 

**** STABILITY DERIVATIVES AS A FUNCTION OF INERTIA (AND PSI) 

SDMU = (HR*(DHDU + THRUST*DADU) + DMDU)/YYI 
SDMQ = (HR*(DHDQ + THRU ST*D ADQ) + DMDQ)/YYI 
SDLV = -SDMU*YYI/XXI 
SDLP = SDMQ*YYI/XXI 

**** STATE MODEL 

UBD = (UB+DEL1)*SDXU + QB*SDXQ - G*THETR 
WBD = (WB +DEL2) * SDZW 
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QBD 

VBD 

PBD 

UB 

WB 

QB 

THETR 

VB 

PB 

PHIR 

UBDP 
WBDP 
QBDP 
QBP 
VBDP 
PBDP 
PBP 


(UB+DEL1)*SDMU + QB*SDMQ + QBDFB 
VB*SDYV + PB*SDYP + G*PHIR 
VB*SDLV + PB*SDLP + PBDFB 


= UB + 0.5*DT2*(3.*UBD - 
= WB + 0.5*DT2*(3.*WBD 
= OB + 0.5*DT2*(3.*QBD - 
= THETR+ 0.5*DT2*(3.*QB 
= VB + 0.5*DT2*(3 *VBD - 
= PB + 0.5*DT2*(3.*PBD - 
= PHIR + 0.5*DT2*(3.*PB - 

= UBD 
= WBD 
= QBD 
= QB 
= VBD 
= PBD 
= PB 


UBDP) 

- WBDP) 

QBDP) 

-QBP) 

VBDP) 

PBDP) 

PBP) 


C 

C 


CO RATES AND POSITIONS RELATIVE TO THE EARTH 

VFD = UB*COS(THETR) + VB*SIN(PHIR)*SIN(THETR) + 
i WB*COS(PHIR)*SIN(THETR) 

- VR*GOS(PHIR) - WB*SIN(PHIR) 

7ED = -UB*SIN(THETR) + VB*SIN(PHIR)*COS(THETR) + 
\ WB*COS(PHIR)*COS (THETR) 

yc - XE + 0 5*DT2*(3.*XED - XEDP) 

W I YE X 0.5*DT2* 3.*YED - YEDP) 
zf = ZE + 0.5*DT2*(3.*ZED - ZEDP) 

XEDP = XED 
YEDP = YED 
ZEDP = ZED 


C 

C 


AERODYNAMIC FORCE AND MOMENT CALCULATIONS, BODY AXIS 


300 


PTX = UBD*XMASS 

PXY = VBD*XMASS 

prpZ = WBD*XMASS 

TTL = PBD*XX1 

TTM = QBD*YYI 

CONTINUE 


RETURN 

FND 


TITLE 
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c 

c 


SUBROUTINE PRTOUT 

SUBROUTINE PRTOUT 


C 

c CREATION/MODIFTCATION LOG: 


3/88 


C 
C 

c 
c 
c 
c 
c 
c 
c 
c 
c 
c 
c 

CIB 

CIPRT 

C 

C 

CHO 

C TIME 

CXTHET 

CXPHI 

CXPSI 

CHITEL 

C HITER 

C 

C 

C 


date name 


REMARKS 
J. TOTAH (NASA) WRITTEN 


INTRODUCTION: 

SSSSgw AIRCRAFT states and stability 
DEHNITIONS OF VARIABLES ARE AS FOLLOWS: 


lAfllj tootflag N prin ™ g every second 

“S S gggig * m0M ™ gr ound 

DA(086) POSITION OF ROTORS IN TPP 
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non 


EQUIVALENCE (DA(009),SDLP ) 

EQUIVALENCE (DA(028),RADIUS ) 

EQUIVALENCE (DA(052),UBD ) 

EQUIVALENCE (DA(053),VBD ) 

EQUIVALENCE (DA(054),WBD ) 

EQUIVALENCE (DA(055),PBD ) 

EQUIVALENCE (DA(056),QBD ) 

EQUIVALENCE (DA(057),UB ) 

EQUIVALENCE (DA(058),VB ) 

EQUIVALENCE (DA(059),WB ) 

EQUIVALENCE (DA(060),PB ) 

EQUIVALENCE (DA(061),QB ) 

EQUIVALENCE (DA(072),THETR ) 

EQUIVALENCE (DA(073),PHIR ) 

EQUIVALENCE (DA(080),XE ) 

EQUIVALENCE (DA(081),YE ) 

EQUIVALENCE (DA(082),ZE ) 

EQUIVALENCE (DA(083),FTX ) 

EQUIVALENCE (DA(084),FTY ) 

EQUIVALENCE (DA(085),FTZ ) 

EQUIVALENCE (DA(086),TTL ) 

EQUIVALENCE (DA(087),TTM ) 

EQUIVALENCE (DA(088),DT2 ) 

EQUIVALENCE (DA(089),PSIR ) 

EQUIVALENCE (DA(100),HITE ) 

EQUIVALENCE (DA(111),AD ) 

EQUIVALENCE (DA(113),CD ) 

EQUIVALENCE (DA(1 17), DELI ) 

EQUIVALENCE (DA(1 18),DEL2 ) 

EQUIVALENCE (DA(123),HO ) 


DATA HO / 10. / 


C**** CODE TO PRINT EACH SECOND 
C 

XPSI = PSIR*57.3 
XPHI = PHIR*57.3 
vtuft = THFTR*57 3 

HITEL = AD+CD+ZE+RADIUS*(PHIR*COS(PSIR)-THETR*SIN(PSIR)) 
HITER = AD+CD+ZE-RADIUS*(PHIR*COS(PSIR)-THETR*SIN(PSIR)) 


IF (IMODE .EQ. -1) IB = 1; TIME = -DT2*2. 

TIME = TIME + DT2 

IPRT = ABS(INT(DT2*1PART) - IB) 

IF (IPRT .EQ. 1 .AND. ICHKDYN ,EQ. 1) GO TO 10 
IF (IPRT .EQ. 0) GO TO 10 
GO TO 20 
10 CONTINUE 
IB = IB + 1 


C 

C**** WRITE STATEMENTS 
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uuu 


WRITE (3,100) TIME 

WRITE (3, 1 10) FTX, UBD, UB, XE , XTHET 
WRITE (3,120) FTY, VBD, VB. YE , XPHI 
WRITE (3,130) FTZ, WBD, WB, ZE , HITE 
(3 ’ 140 ) PBD . PB, DELI , HITEL 

WRITE oS™.’ QBD ' QB ’ DEL2 ’ H ' TCR 
WRITE (3,170) SDXU, SDMU, SDYV, SDLV 
WRITE (3,180) SDXQ, SDMQ, SDYP, SDLP 
WRITE (3,190) SDZW 

IF (ISAS .NE. 0) WRITE (3,200) 

IF (HITER .LT. 0 .OR. HITEL .LT. 0) ID = 1; WRITE (3,210) 
20 CONTINUE 
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190 

200 

210 


FORMAT(////20X, 'AIRCRAFT STATES AT ' F5 2 ' SEC /) 
FOR ™JlZ?' X z!.’ p7 - 2 >' UBD=',F7.2,’’ UB=’,F7!2,' 


VBD=',F7.2,' VB=',F7.2,’ 

WBD=’,F7.2,' WB=’,F7.2,' 


1 ’ XTHET=',F7.2) 

FORMAT('FTY=',F7 2 ' 

1 ' XPHI=',F7.2) 

FORMAT('FTZ=',F7.2,’ 

1 1 HITE=',F7.2) 

FORMAT(TTL=',F7 2 ' 

1 ’ HITEL=',F7.2) 

FORMAT('TTM= , ,F7.2,’ 

1 ' HITER=',F7.2) — » 

5X, 'AIRCRAFT STABILITY DERIVATIVES AT PST - ' F 7 o a 
FORMAT('SDXU<F9.2,' SDMU<F9.2,' SD Yv' F92' SDI V-' F 9 
FORMAT('SDXQ=',F9.2,' SDMO= F9 2 1 SDYP’F) 2 ' mp'coi 
FORMAT('SDZW=’,F9.2) ’ iUYF ’ F9 - 2 ’ SDLP = ,F9.2) 

FORMAT(/20X,’THE CONFIGURATION IS AUGMENTED') 
FORMAT(10X,'******************** CRASH ********************>) 


PBD=',F7.2,' 

QBD=',F7.2,' 


PB=',F7.2,' 

QB=',F7.2,' 


XE=',F7.2, 

YE=',F7.2, 

ZE=',F7.2, 

XDEL1=',F7.2, 

XDEL2=',F7.2, 


RETURN 

END 
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TENSION CABLES 


CENTER SECTION 



ROTOR TIP 


CONTROL SURFACE 



Figure 1 Aircraft description. 
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IMAGINARY IMAGINARY 
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(a) LONGITUDINAL AXIS (a: ^ = 0°, b: ip = 45°. c: \p = 90°) 
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REAL 

(b) LATERAL AXIS (a: ^ = 0°, b : \p = 45°, c: $ = 90°) 

Figure 3.- Unaugmented stability analysis results, a. Longitudinal axis (a: y = 0°, b: y = 45°, 
c: \\f = 90°). b. Lateral axis (a: \|r = 0°, b: \|/ = 45°, c: y = 90°). 
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IMAGINARY 



Figure 4.- Augmented stability analysis results. 
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2(101 + 19l) t= 1 .25-3.75 sec 
1(101 + K9l) t =l 1.25-13.75 sec 





RATE LIMIT, rad/sec 
Figure 5.- Parameter variation results 
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Figure 6.- Continued. 
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TIME, sec 

Figure 7.- Continued. 
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TIME, 


Figure 7.- 
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P, rad/sec 2 V, ft/sec 2 U, ft/sec 




, rad/sec q, rad/sec p, rad/sec W, ft/sec V, ft/sec 



Figure 8.- Continued. 
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, mph u', mph Vc$' rat * 



Figure 8.- 


5 







Figure 10.- Inertial representation. 
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Figure 12.- Discrete simulation flow diagram. 
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